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Abstract 

A simple model is proposed that describes the various morphodynamic principles of migrating 
cells from polar to amoeboidal motions. The model equation is derived using competing internal 
cellular compass variables and symmetries of the system. Fixed points for the N = 2 system 
are closely investigated to clarify how the competition among polaritors explains the observed 
morphodynamics. Response behaviors of cell-to-signal stimuli are also investigated. This model 
will be useful for classifying high-dimensional cell motions and investigating collective cellular 
behaviors. 
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The migration of cells on substrates is a key component of various biological functions, 
such as development, immune system response, and wound healing Cells exhibit diverse 
and involved morphodynamics depending on their type, developmental stage, and environ- 
mental conditions, and yet represent ordered and common dynamics across a range of species 



3j . Many cells exhibit an apparent monopolar shape, consisting of a head and a tail. The 
keratocyte-like crescent shape is also widely observed in nature |4j . The highly investigated 
Dictyosterium cells form a roundish shape with amoeboidal, non-directed, random protru- 
sions. Under poor nutrient conditions, these cells become elongated with a definitive polarity 
and exhibit a zig-zag ( "split-and-choice" ) motion {5, 6], followed by a collective spiral migra- 
tion pattern [7j] . These cellular behaviors are regulated by the cytoskeleton (specifically, actin 
filaments) mediated by signaling molecules (e.g., phosphoinositide lipids). It is now recog- 
nized that interactions among these molecules lead to the instability of a uniform molecular 
distribution inside a cell and gen erate self-organized chemical waves to support complex cell 



morphodynamic processes js-ll3|. A number of theoretical models have been proposed by 
considering associated chemical reactions BOH 

or actin polymerization |19l-l24|. How- 
ever, these models are primarily focused on the onset of instabilities and a relatively simple 
pattern, and not generally intended to elucidate how a variety of morphodynamic processes 
from ordered to amoeboidal cell motions are organized. One of difficulties in addressing this 
issue is the requirement of large computational power for the execution of these models. 
Thus, simplified modeling with appropriate abstraction is another theoretical challenge in 
identifying the mechanism of diversity in cellular morphodynamics. 

Here, we adopt one of the familiar concepts known as a "chemical compass," which was 
introduced as a hypothetical internal cellular state and an intuitive representation of intrinsic 
cell directionality (Fig. H{a)) |27H32|. The compass can be interpreted to represent cellular 
polarity as dictated by molecular distribution. Based on the potential high dimensionality of 
molecular dynamics, there may be multiple cell compasses rather than only a single one. In 
practice, this is evident in amoeboidal motion, which exhibits a number of protrusions with 
patched molecular localization at the cell boundary (Fig. QJb)). This is also supported by 
the observed mechanisms by which a cell responds to signal stimuli from different directions; 
cells turn by rotating their existing head in response to signals, although the existing head 
occasionally disappears in response to signals from the rear and a new head is formed in the 
direction of the signal (Fig. [He)) 0. These observations motivated us to introduce multiple 
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FIG. 1. (a) Polarized cell represented by a single cellular compass variable (polaritor) (b) Amoe- 
boidal cell represented by multiple polaritors. (c) Distinct responses to signal stimuli from different 
directions. The existing cell head turns in response to stimuli from the front to change direction 
(top), whereas the cell head is replaced by a new head in response to signals from the rear (bottom). 

compass variables for describing cellular dynamics. A model is derived in the following 
text for single cellular motion based on the competitive dynamics occurring among the 
hypothetical compass variables. 

Intrinsic compass variables in the cell are represented by complex variables W% = 
Rit x6i (i = 1,2, . . . , N) with their amplitudes (Ri) and directions (9i). We referr to the 
variables as polaritors in this study. A part of polaritors is supposed to be simultaneously 
active (i.e., Ri > 0), while the other polaritors are in a quiescent state (i.e., Ri = 0). 
The dynamics of these polaritors should satisfy the following requirements: First, a sin- 
gle polaritor spontaneously breaks U(l) symmetry, thereby obeying an equation similar to 
W = W — \W\ m W for an even integer m. Second, the system should satisfy all plausible 
symmetries; in other words, the equation should be invariant by the following transforma- 
tion of variables: (i) simultaneous shift of polaritors' direction 9^ — > 9i + 6 (isotropy), (ii) 
reflection 9{ — > —9i for all i (mirror symmetry), and (iii) permutation among polaritors' 
indices. Third, the system should possess invariant subspaces specified by Wi = for each 
i, since polaritors in quiescent states should not influence the dynamics of the system. On 
;he basis of these requirements, the following equation is proposed for the dynamics of Wi 

m 

W i = W i -\W i \ A W i -^ (T,'\Wj\ 2 ) Wi-8 (Z'Wf) Wi 

+ Cl (VWj) \W t \ 2 + c 2 (E'Wj) W? + (1) 

In this equation, Wi is the complex conjugate of Wi, E' = E^ represents the calculation 
of the sum over j except i, and the coefficients 7, 5, c\ and C2 are real numbers due to the 
requirement of symmetry (ii). m = 4 is adopted to ensure the stability of the system for 
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arbitrary parameter values. The final term is added to incorporate white Gaussian 
noise into the statistics (&(£)) = and (t')) = 2DSijS(t — t')/N. In polar coordinates, 

Eq. flTJ, without the noise term, is expressed as 



where p = c\ + c 2 and 1/ = c\ — c 2 . These equations indicate that 7 represents angle- 
independent competition among polaritors, while 5 and p represent angle- dependent inter- 
actions of the first and second order. The velocity of the cell centroid x is set to obey the 
equation 



For visibility, the shape of the cell is assumed to depend on polaritors as an angular-radius 
function L(6) = c^L^Q) with L (6) = (R + J2i Rie Acos{e ~ ei) ) and with the normalization 
factor c 2 = J^ 7T Ll(9)d9 to keep the cell area constant. 

To understand how the competition among polaritors gives rise to new dynamics, we first 
examined the simplest N = 2 system in the absence of noise. With definitions of relative 
angle ip = 9 2 — Q\ and mean angle \l/ = (0i + 9 2 )/2, Eqs. fl2j) and (j3J) reduce to 



The equations for Ri,R 2} and ip are incorporated within these variables; ^ evolves depending 
on these variables, and thus, the dynamics related to (R\, R 2 ,ip) is the primary interest 
herein. The domain of ip can De restricted to — tt < ip < tt. We observed that a variety 
of fixed points appear in this system (Fig. E^a)) as listed below, where several obvious 
fixed points obtained by replacing % — 1 and 2 are excluded. The stability of these fixed 
points are also determined by considering the linear equation p = Jp, where p indicates a 
small deviation in (Ri, R 2 ,ijj) and J is a Jacobian matrix at a fixed point. Eigenvalues and 



Ri = Ri- Rf - £' (7 + S cos 2(0,-0*)) R)R { 

+ pY! co^{Qj-e,)RjR\, 
9\ = -8Y.'R) sin 2(0,-0*) + vTt sin (0,-0,)^^ 



(2) 
(3) 




(4) 



Ri = Ri -Rf - (7 + 5 cos 2tp)R 2 2 R l + p cos ipR 2 R\, 
# 2 = # 2 _#5_( 7 4. Scos2tfj)RjR 2 + p cos tfjRxRl, 

ij; = 5(Rl + R 2 2 ) sin 2i\) - 2v sin tpRiR 2 , 

if = 5{Rl-Rl) sin 2^/2. 



(5) 
(6) 
(7) 
(8) 
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FIG. 2. (a) Fixed points of the N = 2 system and their relationship (see text). Red and blue 
arrows indicate polaritors. (b) Phase diagram of the N = 2 system. 7 and c\ are fixed as 
7 = 0.7 and c\ = 1.0. Different types of fixed points appear, depending on the initial condition 
in regions indicated by * (Ao and So) and ** (So and S^). P n may appear when 7 is sufficiently 
large (7 > 1 + \5\). (c-f) Examples of numerical simulations for the N = 2 system. Red lines 
indicate the trajectories of cell centroids. Parameters (7, 5, 01,02) are set as (c) (1.2,0.1,1.0,0.5), 
(d) (0.7,0.4,1.0,1.0), (e) (0.7,0.2,1.0,1.6), and (f) (0.78,0.38,-0.45,0.45). Parameters for shape 
are set as i?o = 1-8, A = 2.4, and A = 0.3. 



corresponding eigenvectors of J are denoted as A = (Ai,A 2 ,A 2 ) and pi,p 2 ,p*3, respectively. 
Fig. |2(a) summarizes the fixed points and their relationships. 

[Z] (Ri,R 2 ) = (0,0) for arbitrary ip. This is a trivial fixed point with Ai = A 2 = 1 and 
A3 = 0, and always remains unstable. 

[P n ] (Rx,R2,ip) = (1,0, twt/2) with n as an integer (Fig. 12(a)). These fixed points are 
also an expected trivial state, where one polaritor is active while the other is quiescent. 
Eigenvalues are denoted as Ai = — 4, A 2 = 1 — (7 + 8cos2%l)) ) and A 3 = 25 cos 2tp. One of 
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these fixed points is stable as long as 7 > 1 + 

[So, S^, S s ] Fixed points with symmetric finite amplitudes (Ri = R 2 = R>0; Fig. [2(a)). 
Directions of the two polaritors coincide for So (ip = 0), but are reversed for S n (ip = 7r). S s has 
a non-trivial relative angle ip s that is determined by cos^ = v/25. Square amplitudes Ro, 
R w , and R s are defined by R 2 = [—(7 + 5 cos 2ip — /i cos ■?/>) + a/(7 + 5 cos 20 — [icosvp) 2 + 4]/2 
with the substitution of ip — 0, 7r, and T/> a , respectively. The linear stability of these fixed 
points are evaluated by Ai = 2 — 6i? 4 , A2 = — 2 — 2i? 4 , and A3 = 2R{25 cos 20 — z/ cos ■0). 
The eigenvector of A 3 is p 3 = (0,0,1), indicating that A 3 defines stability in the direction 
of the relative angle ip (angular stability). So and S^ are unstable in the angular direction 
in the regions v < 25 and v > —25, respectively, in which S s appears through pitch-fork 
bifurcation. In contrast, Ai has the eigenvector p\ = (1,-1,0) and defines stability in the 
amplitudal directions; Ai > indicates a break in the amplitudal symmetry. A2 is always 
negative and does not alter stability. 

[Ao] Fixed points at which two polaritors have distinct finite amplitudes and the same 

direction (Ri> R 2 >0 and ip = 0; Fig. [2(a)). The condition of the fixed points is defined by 

R\+R 2 =^/+5 and Rf+R 2 +('-f + 5) 2 = 2fiR 1 R 2 (see Eqs. (JMZ]))- Two types of solutions appear 

as a pair of saddle-node bifurcations at 4(7 + 5) 2 + fi 2 = 4. One of them is always unstable, 

irrespective of parameter values, and is never realized. Only the other fixed point, denoted as 

A , can be realized. By denoting 7' = 7+5, square amplitudes of the solution are defined as 
i?2 = [y + {y2_(_ /i _ v / At2+47 ,2_ 4 )2|i/2 ]/2 and R 2 = [y _{y2 _^2 +4 y2 _ 4)2}l/2]/ 2 . 

Through the change in the amplitudal direction, this fixed point connects to So via pitch-fork 
bifurcation and connects to P n via transcritical bifurcation (Fig. |2]^a)). Angular stability is 
lost when the eigenvalue A3 = 27'^ — 2vR\R 2 is positive, with an eigenvector of p*3 = (0,0, 1). 

[A,r] Fixed points at which two polaritors have distinct, finite amplitudes and opposite 
directions [R\ > R 2 > and ip = n; Fig. [2]^a)). Similar to the case of A , two types of fixed 
points appear as a pair of saddle-node bifurcations at 47 /2 + ji 2 = 4. One of these points 
is always unstable, and only the other, A,,-, can be realized. The solution relates S^ and 
P n via pitch- fork and transcritical bifurcations, respectively (Fig. [2(a)). Angular stability 
is lost when A3 = 25(7 +5) + 2vR\R 2 is positive. Square amplitudes at A^ are defined as 
R\ = Yi' +{ 7 '2 -( /U - v / / u2 + 4 7 ' 2 -4) 2 } 1 /2]/ 2 and R 2 = [i -{ 7 ' 2 -(^ - v ^ 2 +ly 2 ^4) 2 } 1 / 2 ]/2. 

[A a , Ab] Fixed points with distinct, finite amplitudes (Ri > R 2 > 0) and a ip that is 
neither nor tt (Fig. [2(a). These solutions are obtained from the following conditions derived 
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from Eqs. flM2): R\ + R\ = 7 + <5cos2^, i2f + R\ + (7 + <5cos2^) 2 = 2{icos^R 1 R 2} and 
5 {R\ + R\) cosip = vR\R 2 . The expression giving the solutions is lengthy and of secondary 
importance for the purpose of this work. We showed that there can be twelve types of 
possible solutions. Two of them bifurcate from Ao and A,,- by the angular instabilities 
through pitch-fork bifurcations. Let us denote them as A a and Ab, respectively. These 
two fixed points are separated, since Ri = R 2 = at ip = ±7t/2; is less than tt/2 for A a 
but larger than tt/2 for Ab- The other fixed points are mostly unstable; two of them may 
appear in the limited parameter regions via saddle-node or subcritical bifurcation (these 
fixed points are not discussed further in this report). 

Figure EJ^b) shows the phase diagram against c 2 = (/U — v)/2 and 5, with fixed values of 
7 = 0.7 and c 2 = 1.0. At the fixed points P n , S , S s , A , and A n , cells move directionally 
straight = 0), but they have different shapes. Examples of cell motion at P n and S s are 
shown in Fig. |2]^c) and (d), which are similar to polarized and keratocyte-like cell motions, 
respectively. On the other hand, \1/ is nonzero and the cells show migration with a circular 
orbit at the fixed points A a and A b (Fig. EJ^e)). A cell at fixed point S n is bipolar and does 
not exhibit migration. 

In addition to these fixed points, numerical simulations revealed oscillatory dynamics 



in narrow parameter regions |34j. The cell in the oscillatory state exhibits a quasi-periodic 
Lissajous orbit, as shown in Fig.[2](f). Another type of motion found by numerical simulation 
is repetitive right angle turns (Fig. E^g)). Similar types of motions were reported for another 



model 



26). 



Taken together, many inner states appear on account of the competition among polaritors, 
even in the simplest system of N = 2. Note that Eq. ([7]) indicates that ip = and ip — 7T 
are separatrices in the phase space of the N = 2 system. In addition, r = R 2 — R 2 obeys 
the equation of the form dr/dt = Q x r from Eqs. (JSJ) and fl6]), where Q is a function of 
(Ri, R 2 ,ip), indicating that R\ = R 2 is an additional separatrix in the system. Thus, the 
order in amplitudes of polaritors defined by the initial condition is maintained in a noiseless 
system. These separatrices constrain the dynamics of the N = 2 system; for example, zig-zag 
motion is forbidden. 

Such constraints are absent for N > 2; therefore, the system can exhibit various types 
of motion. For N = 3, zig-zag motion of cell migration is observed, as shown in Fig. Efa), 
where two oscillating polaritors dictate the position of the head of the cell and periodically 
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FIG. 3. (a) Zig-zag motion for the N = 3 system with parameters (7, 5, c\, 02) = (0.78, 0.3, 0.6, 0.8). 
(b) Chaotic motion for the N = 3 system. Parameters are set as (0.8,0.55,0.35,-0.8). (c,d) 
Examples of chaotic motions in the N = 20 systems. Parameters are set as (0.7, 0.3, 0.6, 0.8) (c) 
and (1.2, -0.22, 0.8, 2.35) (d). Parameters for shapes are set as (Rq,A, A) = (1.5, 12.0, 0.3) for (a,b) 
and (1.8,4.0,0.3) for (c,d). 



determine the direction of migration. Chaotic inner dynamics is found in the case where 
the cellular trajectory is also a chaotic orbit, as shown in Fig. El^b). For larger N systems, 
additional types of dynamics appear. Fig. [3](c) shows an example of cellular motion at the 
N = 20 system, wherein the cell shape can fluctuate significantly, similar to the case of 
amoebic motion. Fig. [3](d) shows another example in the N = 20 system, for which the cell 
exhibits repetitively straight motion, followed by locally diffusive random migration. 



Finally, we considered the chemotactic behavior of a cell in a signal field S(pc, t) 35|, 136 ]. 
Cell shapes and internal cellular compasses are correlated with the sensing ability to external 
signal molecules. Assuming that a cell can sense the gradient of the signal field s(x) = 
VS*(x) = s(x)e"^ x \ the coupling of polaritor variables with the signal field is incorporated 
by combining dis\ Wi\ 2 + d^Wf with Eq. ([T|). For a = di + 62 and k = d\ — d 2 , the following 
additional terms appear in Eqs. (T5]) and (EJ): 



sig 



asR 4 cos(8j — 



QSig 



(9) 
(10) 



nsRi sin(0 — 9i). 

Here, both a and k are set as positive for an attractant signal. For a single polaritor with 
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FIG. 4. (a) Set-up of the numerical simulation. A cell with a polaritor directed in the x-direction 
receives a signal from direction <p at t = 0. (b) Probability of cell directional change by the 
replacement of polaritors is shown in a color scale against the invoked signal direction if) and 
parameter k. The replacement can occur only when <j) > ir/2. Parameters are set as 7 = 1.4, 
S = 0.15, n = 0.3, v = —0.1, a = 2.0, and s = 1.0. A small amount of noise is added during the 
simulation (D = 5.0 x 10~ 3 ). The probability is calculated from 100 independent simulations for 
each parameter set. 

J?~l, the direction develops as 9{ ~ sin ((/> — #;) and is directed to the maximum gradient 
of the signal concentration, as introduced previously [37]. When \6i — <j>\ is larger than 7r/2, 
Rf 9 becomes negative and the existing polaritor begins to shrink. Thus, the replacement of 
polaritors may occur, depending on the direction of the signal gradient. 

To confirm the response of cellular behaviors, numerical simulation is conducted with 
the following settings. First, a cell with N = 2 polaritors at the single polarity state P n 
is prepared in the absence of a signal field. The active polaritor is set in the x-direction. 
Then, a signal field is applied at t = with s(t) = se^Qit), where s and are constants and 
0(-) is a Heaviside function (Fig. 11(a)). Weak noise is added (D = 5.0 x 10~ 3 ). Then, the 
system shows either rotation of the existing polaritor (Fig. H(c), top) or the replacement of 
polaritors (bottom), depending on signal direction <ft and other parameters. Fig. 11(b) shows 
the probability of the replacement occurring for respective signal direction tp and parameter 
k. The results demonstrate that the model cell can respond to signal stimuli from the rear 
(<f) > 7r/2) by switching polaritors. 

By considering competition among polaritors, the present study investigates a mecha- 
nism for organizing a variety of cellular behaviors linked to morphology and migration. The 
proposed model exhibits distinct polar, keratocyte-like, zig-zag, and chaotic amoeboid mo- 
tions that are relevant to experimental observations. This model is only constrained by the 
symmetry of the system. As demonstrated in different systems like quadrupedal locomotion 
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such an approach based on symmetry can be helpful in classifying a variety of possible 
cellular motions. In addition, the model is quite simple and requires little com put ational 
power, making it possible to use the model to study collective cellular behaviors [39|, |40 |. 

An advantage of our model is that it is easily extensible to higher dimensions even with 
the same number of model parameters, and it provides an intuitive and consistent interpre- 



tation of cellular behavior. Previous models have been reported |25|, |26J that exhibit similar 
dynamics to those presented in this study. However, these models become complicated by 
including higher-order tensor variables; in fact, zig-zag and chaotic motions appear in the 
equations that contain more than 20 parameters [26]. 

Because the concept of a polaritor is introduced here as a rather abstract variable, a 
future step will be to identify the molecular basis of the polaritors and their interactions. 



Va 
6. 



idating Eq. (JT]) from detailed subcellular processes (e.g., reduction from detailed models 
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24|) will elucidate the way in which cellular motion depends on molecular parameters, 



which improves the correspondence of the model with experimental observations. 
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